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1. Introduction 

In this paper we discuss how primordial non-Gaussianity (NG) in the cosmological 
perturbations is left imprinted in the Large-Scale Structure (LSS) of the universe in a 
ACDM cosmology. We show how the information on the primordial non-Gaussianity, set 
on super-Hubble scales, flows into smaller scales through a complete General Relativistic 
(GR) computation. Primordial NG thus leaves an observable imprint in the LSS. 
Another interesting finding is that, on sufficiently large scales, there is another additional 
source of non-Gaussianity which arises from GR corrections, the leading contributions 
of which are Post-Newtonian terms as first pointed out in Ref. pQ. The importance 
of the signatures of primordial non-Gaussianity in the evolution of the matter density 
perturbations is due to the fact that future high-precision measurements of the statistics 
of the dark matter density will allow to pin down the primordial non-Gaussianity, thus 
representing a tool complementary to studies of the Cosmic Microwave Background 
(CMB) anisotropies. It is beyond the scope of this paper to go into the details of the 
various theoretical and observational methods related to this issue. This paper will 
then serve as a basic guideline to capture the starting point expressions that relates 
primordial NG to the dark matter density and gravitational potentials, outlining some 
specific non-Gaussian signatures in the LSS that can be potentially interesting. We 
perform our calculations in assuming a flat Universe with pressure-less matter, i.e. Cold 
Dark Matter (CDM) plus non-relativistic ordinary matter and a cosmological constant, 
hereafter ACDM cosmology. We present our results both in the Poisson gauge and in 
the comoving time-orthogonal gauge, which are relevant for comparison to observations. 

The primordial NG considered in our analysis is set at primordial (inflationary) 
epochs on large (super-Hubble) scales. At later times cosmological perturbations reenter 
the Hubble radius during the radiation or during the matter- and dark energy-dominated 
epochs. For scales re-entering the Hubble radius during the radiation dominated era one 
should include the radiation in the evolution equations, thus using a complete second- 
order matter transfer function also for those scales. A detailed treatment of it has 
been given in Refs.[2j [3]. Refs. jU |5] also investigate the evolution of the dark matter 
perturbations up to second-order accounting for a radiation-dominated epoch, both 
analytically and numerically. Here we will focus on large scales for which the effects 
arising during the radiation-dominated epoch can be neglected. 

The plan of the paper is as follows. In Section 2 we give the general form of the 
perturbed line-element and we introduce the gauge-invariant curvature perturbation 
of uniform density hypersurface at first and second order, which is used to provide 
the inflationary initial conditions, including the effect of primordial non-Gaussianity. In 
Section 3 we derive the second-order expression for the density perturbation in a ACDM 
cosmology, in the Poisson gauge, taking into full account both NG initial conditions 
and Post-Newtonian corrections arising from the non-linear evolution of perturbations 
according to the fully General Relativistic treatment. Section 4 contains a similar 
calculation in the comoving-synchronous gauge. Section 5 contains our concluding 
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remarks. 

2. Metric perturbations and primordial non-Gaussianity 

We consider a spatially flat Universe filled with a cosmological constant A and a non- 
relativistic pressureless fluid of Cold Dark Matter (CDM), whose energy- momentum 
tensor reads T^ v = pu^u u . Following the notations of Ref. [6], the perturbed line element 
around a spatially flat FRW background reads 

ds 2 = a 2 (r){-(l + 2<p)dr 2 + 2u i drdx i + [(1 - 2^)<% + XijWd^} . (1) 

where a(r) is the scale factor as a function of conformal time r. Here each perturbation 
quantity can be expanded into a first-order (linear) part and a second-order contribution, 
as for example, the gravitational potential <p — + (fi^/2. Up to now we have not 
choosen any particular gauge. We can employ the standard split of the perturbations 
into the so-called scalar, vector and tensor parts, according to their transformation 
properties with respect to the 3-dimensional space with metric Sij, where scalar parts 
are related to a scalar potential, vector parts to transverse (divergence-free) vectors and 
tensor parts to transverse trace-free tensors. Thus <ft and ip, the gravitational potentials, 
are scalar perturbations, and for instance, = diU^ + uj> t , where wW is the scalar 
part and is a transverse vector, i.e. d^uf = ((r) = (1,2) stand for the rth- 
order of the perturbations). The symmetric traceless tensor Xij generally contains a 
scalar, a vector and a tensor contribution, namely Xij = ^ijX + d%Xj + djXi + Xij: 
where = didj — (l/3)V 2 5ij, Xi is a solenoidal vector (d t Xi = 0) and x%j represents 
a traceless and transverse (i.e. d l x%j = 0) tensor moddjl As for the matter component 
we split the mass density into a homogeneous p(r) and a perturbed part as p(x, r) = 
p(r)(l +S^' + 5^/2) and we write the four velocity as = (5$ +t> /i )/a with u^u^ = — 1 
and = v^ 1 ^ + 

The Friedmann background equations are 

3^ 2 = a 2 (87rGp(r) + A), (2) 

and 

p'{r) = -m P {r) , (3) 

where a prime stands for differentiation with respect to conformal time, and T-L = a' /a. 
The matter density parameter is fi m (r) = 87rGa 2 (r)p(r)/(3'H 2 (T)). 

Before recalling how one can parametrize the primordial non-Gaussianity, we need 
to provide a general definition for the amplitude of non-Gaussianity characterizing the 
matter density contrast beyond the usual second-order Newtonian contributions. It 

J In what follows, for our purposes we will neglect linear vector modes since they are not produced 
in standard mechanisms for the generation of cosmological perturbations (as inflation), and we also 
neglect tensor modes at linear order, since they give a negligible contribution to LSS formation. 
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proves convenient to introduce an effective gravitational potential obeying the Poisson 
equation 

- v 2 $ = ^n m n 2 s, (4) 

and at an initial epoch, deep in matter domination, we write 

$in = C+/NL(C 2 -(^ )2 )), (5) 

with the dimensionless non-linearity parameter /nl setting the level of quadratic non- 
Gaussianity, and $ oc g(r), g{r) being the usual growth suppression factor (see Sec. [3]). 
Therefore we will write the matter density contrast in Fourier space in terms of the 
linear density contrast by defining the kernel /Cj(ki, k2; r) depending on the wavevector 
of the perturbation modes as 

4(r) = tf ) (r) + ^ ) W (6) 
We can write the kernel as 

3 g in k 2 

K s (ki, k 2 ; r) = /Cf (ki, k 2 ; r) + -fi m ^ 2 / N L(ki, k 2 , r) m 2 2 , (7) 

2 p(r) 

where k 2 = |ki + k 2 | 2 and /C^ (k 1; k 2 ; r) is the second-order Newtonian kernel. 



2.1. Primordial non-Gaussianity and initial conditions 

We conveniently fix the initial conditions at the time when the cosmological 
perturbations relevant for LSS are outside the horizon. A standard and convenient 
way to account for any initial primordial non-Gaussianity is to consider the curvature 
perturbation of uniform density hypersurfaces £ = Q 1 ' + C,^ /2 + • • •> where = 
— i/j^ — HSp^ 1 ' / p' and at second-order [121 E] 



+ J " «' - ™?) , (8) 

where = ^ + V 2 x (r) / 6 - Tnis 

is a gauge-invariant quantity which remains constant 
on super-horizon scales after it has been generated during a primordial epoch (and 
possible isocurvature perturbations are no longer present). Therefore, provides 
all the necessary information about the primordial level of non-Gaussianity. The 
conserved value of the curvature perturbation ( allows to set the initial conditions for the 
metric and matter perturbations accounting for the primordial contributions. Different 
scenarios for the generation are characterized by different values of (^ 2 \ while the post- 
inflationary nonlinear evolution due to gravity is common to all of them [131 El HS1 Ej- 
For example, in standard single-field inflation is generated during inflation and its 
value is C (2) = 2 (C (1) ) + O (e,rj) pH EH Q3], where e and r\ are the usual slow-roll 
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parameters. Therefore it turns out to be convenient to parametrize the primordial 
non-Gaussianity level in terms of the conserved curvature perturbation as in Ref. [20] 

C^ = 2a NL (C (1) ) 2 , (9) 

where the parameter onl depends on the physics of a given scenario, and in full generality 
it can depend on scale and configuration. For example in the standard scenario onl — lj 
while in the curvaton case (see e.g. Ref. [7]) onl = (3/4r) — r/2, where r « (p a /p)i) 
is the relative curvaton contribution to the total energy density at curvaton decay. In 
the minimal picture for the inhomogeneous reheating scenario, onl = 1/4. For other 
scenarios we refer the reader to Ref. [7J [8j [9j HOI [EE]. One of the best techniques to 
detect or constrain the primordial large-scale non-Gaussianity is through the analysis of 
the CMB anisotropies, for example by studying the CMB bispectrum [3, [181 [IS] • The 
non-linearity parameter /nl as defined in Eq. (jSJ) is defined also to make contact with 
the primordial non-Gaussainity entering in the CMB anisotropies. For large primordial 
non-Gaussianity, when |cinl| ^> 1) /nl — 5qnl/3 (see Refs. [7| 120]). 



3. Dark matter density perturbations at second-order: Poisson gauge 

The goal of this section is to compute the matter density contrast in the Poisson 
gauge [21], i.e. the generalization beyond linear order of the longitudinal gauge, by which 
a more direct comparison with the standard Newtonian approximation adopted in the 
interpretation of LSS observations and in N-body simulations in Eulerian coordinates 
s possible (see, however, Ref. [22] for a critical discussion of the potential problems 
connected to the use of the Poisson gauge on scales comparable with the Hubble radius) . 
In the Poisson gauge one scalar degree of freedom is eliminated from the got component 
of the metric, and one scalar and two vector degrees of freedom are eliminated from g^. 

Let us briefly recall the results for the linear perturbations in the case of a non- 
vanishing cosmological A term. At linear order the traceless part of the (i-j)-components 
of Einstein equations gives <^> = tp^ = if. Its trace gives the evolution equation for 
the linear scalar potential ip 

cp" + m<p' + a 2 Aip = . (10) 

Selecting only the growing mode solution one can write 

<p(x,T)=g(T)(p (x), (11) 

where <po is the peculiar gravitational potential linearly extrapolated to the present time 
(to) and g{r) = D + {j)/a{r) is the so called growth-suppression factor, where D + {r) is 
the usual linear growing-mode of density fluctuations in the Newtonian limit, i.e. the 
no-decaying solution of the differential equation 

D" + HD' — ^-T-i 2 Vt m D = . (12) 
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We normalize the growth factor so that D + (tq) = a$ = 1. The exact form of g can be 
found in Refs. [231 12H 125]. In the A = case g — 1. A very good approximation for g 
function of redshift z is given in Refs. [23l El] 

-l 



with fL 



fWi + ^) 3 /£ 2 ' 



(l + fi m /2)(l + fi A /70)J , (13) 

; mAi -r (zj, = n 0A /E 2 (z), E{z) = (1 + z)U{z)/U Q = 

[Q 0m (l + z) 3 + ^oa] 1 ^ 2 and f2 0m , fi A = 1 — ^om, the present-day density parameters 
of non-relativistic matter and cosmo logical constant, respectively. According to our 
normalization, g(z = 0) = 1. The energy and momentum constraints provide the 
density and velocity fluctuations in terms of ip (see, for example, Ref. [13] and [26l 
for the A case) 

1 



SO) 



X 1 ) 



ATiGa?p 
1 



- [VV - (if' + Hep)] 



(14) 

^Ga?p^ yT ' ■ (15) 

In a similar way the expression of the second-order matter density contrast can be 
computed starting from the energy constraint given by the (0 — 0) Einstein equation, 
once the evolution of the gravitational potentials is known. The latter has been already 
computed in detail in Ref. [28J and in the following we summarize the main results. 

The evolution equation for the second-order gravitational potential tp^ is obtained 
from the trace of the (i-j)- Einstein equations [§| 

^(2)" + 3^(2)' + a 2 A4j (2) 



S(r) 



(16) 



where S(t) is the source term 
S(r) = g 2 a, n H 2 



if 



(17) 



where, for simplicity of notation, we have introduced 



and 



a {xj 



/(a 



V" 2 (tftpodiifo) - (d'epodjepo) 



dlnD, 



1 



9\r 



d\na ' Hg{f) ' 

which can be written function of Q m as /(fi m ) f 
V -2 stands for the inverse of the Laplacian operator. 



^m(^) ^ 



(19) 
23]. In Eq. ([TJ 



§ The second-order perturbations of the Einstein tensor G M „ can be found for any gauge in Appendix 
A of Refs. [HI [7] and directly in the Poisson gauge, e.g., in Refs. |3]. The perturbations of the 
energy-momentum tensor up to second order in the Poisson gauge have been computed in Ref. |13j for 
a general perfect fluid (see Ref. [7] for expressions in any gauge). 
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The solution of Eq. (1161) is then obtained using Green's method with growing and 
decaying solutions of the homogeneous equation ip + (r) = g(r) and ip-{r) = H(t) / a 2 (r) , 
respectively. The second-order gravitational potentials then read [28] 



^ (2) (t) = [ B x {t) - 2g(r)g in - ^f(a NL - l)g(r)g- m } </?o + ( B 2 {r) - ^g(r)g in ) a (x 



10 



(20) 



(l 



10 



Si(r) + V(r) - 2g(r)g in - — (a NL - l)^(r)^ in ^ + 



5 2 (r) + \g 2 {r) 



x 



e r 



« (x) + E 3 (r)V- 2 a^'(aVo5^o) + S 4 (r)aVo9^o 



(21) 

where we have introduced -Bj(r) = "Hq 2 (/o + 3Qo m /2)~ 1 Bi{r) with the following 
definitions 



Bi(r) 
B 2 (r) 

Bs(r) 
and 

with 



dfU\f){f{f)-\yC{r^) 



dfH 2 (?) 



dr ( e(f 



2(/(r) - I] 



C(r, f) 



■ 3fi m (r) 
5 4 (r) = 



C(r, f ) 



C(r,f) 



e(fi m ) = 



/(f)a(f) 



/ 2 (On) 



g(r)H(f)-g(f 



a 2 (f) 



a 2 (r) 



dfC(r, f) 



H(t) 



(22) 
(23) 
(24) 

(25) 
(26) 



The expression for (j)^ is obtained from the relation between ip^ and (j)^ 



V 2 V 2^(2) = V 2 V ^^_ 4 ^ V 2 V ,^_^2^ + -J ^V 2 (9^ 9Vo)-39^(9Vo9^o)j ,(27) 

which follows from the traceless part of the (z-j)- component of Einstein equations [28J. 
Here (p- m = gi n (fo represents the initial condition taken at some time r in deep in the 
matter dominated era on super-horizon scales. The solutions (J20l and (J2~TT) properly 
account for the non-Gaussian initial conditions parametrized by (onl — !)• These are 
obtained using the expression for during the matter-dominated epoch together with 
Eq. ( l2"7|) and the second-order (O-O)-component of Einstein equations (both evaluated 
for a matter-dominated epoch), so that one can express 0^ and tjj^ in terms of of 
Eq. 0, where C (1) = -5^ in /3 (see Refs. pH ESI d EDI EH]). 

Before proceeding further notice that in the expression for the second-order 
gravitational potentials of Eq. ( 120]) and ( l2Tj) we recognize two contributions. The term 
which dominates on small scales, [^(r) V~ 2 <9j<9 J (<9Vo<9j<^o) + -B^r^VoC^o], which 
gives rise to the second-order Newtonian piece and is insensitive to any non-Gaussianity 
in the initial conditions. The remaining pieces in Eqs. (1201) and (T2TT) correspond to 



2v72a(2). 



4 
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contributions which tend to dominate on large scales with respect to those characterizing 
the Newtonian contribution, and whose origin is purely relativistic. In particular these 
are the pieces carrying the information on primordial non-Gaussianity. For a flat matter- 
dominated (Einstein-de Sitter) universe g{r) = 1 and B\(t) = B 2 {r) = 0, while 
B%{t) —> (5/21)r 2 and B^t) — > — r 2 /14, so that one recovers the expressions of [3] (see 
also [29]). In Ref. [27] second-order cosmological perturbations have been computed in 
the A ^ case from the synchronous to the Poisson gauge, thus extending the analysis 
of Ref. [6], and the CMB temperature anisotropies induced by metric perturbations 
have been also considered by applying the expressions of Ref. [30] (see also Ref. [28]). 
However, an important point to notice is that both Refs. [HI [30] and Ref. [27] disregard 
any primordial non-linear contribution from inflationjjj] 

The matter density contrast at second-order can now be calculated from the (0 — 0) 
Einstein equation 

m^ 2)l + m 2 ^ - vV 2) + i2^ (1 V (1) ' - 2V 2 (^ (1) ) 2 - uu 2 (</> (1) ) 2 - i2u<p w ^ iy 

(28) 



-3U 2 Q r 



We thus arrive at the density contrast in the Poisson gauge 



5® = — 



(/ 



2 A'(t) 



4 



:99b 



+ 



3tt m g 2 
4/V 

3 

1 

V 2 
B 4 (r 
H 



3 

BAi 



g 2 n 

3 



g 2 



l 



Q m g 



B' 2 (r) 
H 



Ag' 



3n 9in 



a(x, T) 



299h 



7^( a NL 



2 Bjjr) 
V m g 2 H 

- l)99m + 2# 2 



vy 

U 2 
- 3B 2 (t 
Bs(r) 



g 2 - m 2 B^r) 



V 2 («9V^), 



u 2 

8 cpV 2 (p 



3tt m g 2 
2 



3Q*m9 2 

^gg-m - 



^m,g 2 

B2(r) 



n 



diipd l tp 



^gg-m 



3Q m g 2 H z 



d i df{&ipd j ip) 



3H 2 B 3 (r) 
2 



3Q m g 2 



(29) 



where A(t) = [-£>i(t) —2gg in — (10/3)(onl — l^^in] an d a ( x ) T ) has the same expressions 
as ao(x) introduced in Eq. (TT8l) with tp in place of <^o- Notice that we have explicitly 
verified that Eq. ( 1291) in the limit of an Einstein-de Sitter universe recovers the expression 
for the matter density contrast obtained in Ref. pQ. 

Eq. (|29|) is the main result of this section. It shows how the primordial NG, which 
is initially generated on large scales, is transferred to the density contrast on subhorizon 
scales. The expression for the density contrast is made of three contributions: a second- 
order Newtonian piece (the last two terms in Eq. (|2"91 . proportional to r 4 in an Einstein- 
de Sitter universe) which is insensitive to the non-linearities in the initial conditions; a 
Post-Newtonian (PN) piece (related to two gradient terms of the gravitational potential) 
which carries the most relevant information on primordial NG; the super-horizon terms 



The results in Refs. [6l [30l [27] have initial conditions corresponding to our onl = 0. 
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(corresponding to the first two lines of Eq. (1291) ). Our findings show in a clear way that 
the information on the primordial NG set on super-Hubble scales flows into the PN 
terms, leaving an observable imprint in the LSS. Another interesting result which shows 
up in Eq. (!29|) is the presence in the Post-Newtonian term of contributions different 
from the primordial NG which are due to weakly non-linear corrections which switch on 
when the modes cross inside the Hubble radius and which can constitute an interesting 
additional source of non-Gaussianity since they can probe large-scale GR corrections. 



3.1. The non-linearity parameter /nl 



We can rewrite the matter density contrast in Fourier space in terms of the linear density 
contrast as in Eq. ([6]). We find 

(k.k 1 )(k-k 2 ) , „ ,_ N A; 2 (ki.k 2 " 



/C?(ki,k 2 ;r) 



3 O 

4 g 2 



Bs(r) 



and for the non-linearity parameter 

|(a NL -l) + l-^- 
3 gi n 



J NL 



(ki,k 2 ; 



+ 



(a NL - 1) + ] 
k 2 ) [2 

(k-k 1 )(k-k 2 



9_ 

Pin 

3 



k 4 

(k-k 1 )(k-k 2 ) k\ 



-H 2 

2 gg\ 



1%) 
2 gg-m 

1 p 2 

— ^_ _l — 

2 Pin 3 



h 2 h 2 
(k 2 



BAt 



h 2 h 2 



,(30) 



kl) 



1 B 2 (t) 

2 gg-m 

3 ^ 2 (r) 
2 



Pin 
2 PPin 



(31) 



ppi 



/lo 3 



k 2 k\k\ 2 PPin ^1^2 

The non-linearity parameter /nl is defined via Eq. (j3J) and in this way it generalizes the 
standard definition of Ref. [31] inferred from the Newtonian gravitational potential. In 
order to obtain the expressions in Eqs. (!30l and (!3~T|) we have performed an expansion in 
(H/hip) < 1 up to terms (H/h) 2 starting from Eq. (125)1 . In Eq. (J3TJ) the primordial NG 
is clearly evident in the piece proportional to (onl — !)• The remaining terms are due to 
the horizon scale Post-Newtonian corrections we commented about in Eq. (129]) . The non- 
linearity (NG) induced by these terms show a specific and non-trivial shape dependence 
that can help in detecting them. Of course their relative importance increases with the 
scale, and in fact it has been shown in Ref. [32] that, through the large-scale halo bias 
techniques investigated in Refs. [331 [M], these GR corrections are potentially detectable. 

It is also worth noticing that some of the terms entering in the non-linearity 
parameter vanish in the limit of a vanishing cosmological constant (e.g. Bi(r) and 
B 2 {t) go to zero in this limit). 

A final comment on Eq. (JBl). In this expression, when the various modes are 
well inside the horizon, one can take the usual configuration for the linear density 



BAt) 



(ki ■ k 2 



~~\~ k<2 3 



BAt) 



2 gg 



perturbations with 5^ (r) oc (2A; 2 T(/c)/3fi m0 ifg)_D + (a), where T{k) is the usual linear 



matter transfer function. In this way one partially accounts for the effects of the 
transition from a radiation- to a matter-dominated epoch. In fact a full computation 
would require to study up to second-order the evolution of the density perturbations 
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also during the radiation dominated era. In this way one can recover a full matter 
transfer function up to second-order. Details about such a computation can be found 
in Refs. [31 IH E]. For example, Ref. [5] shows that the corrections from the full matter 
transfer function, when accounting for a matching at second-order to the radiation 
epoch, gives a relative correction with respect to the Newtonian kernel that is of order 
a eq ~ 10 -4 and that such a correction is equivalent to the effect of a primordial NG 
of /nl ~ 4. That these small-scale corrections are tiny is easy to understand. By 
looking at Eq. ( 128]) one realizes that accounting for the radiation epoch at second-order 
the matter density perturbations from, say approximately the matter-radiation equality 
epoch onwards, will get a correction which scales like 5 ( - 2 - ) | q U i v r 2 , which rescales the 
matching initial conditions to the radiation epoch. However we expect such term to be 
negligible w.r.t. the Newtonian part which scales like {S^\t)) 2 oc r 4 , and also w.r.t. a 
sizable primordial NG. 

4. Dark matter density perturbations at second-order: 
comoving-synchronous gauge 

Let us now see how the second-order matter perturbations and gravitational potentials 
are obtained in the comoving-synchronous gauge. 

We make use of the formalism developed in Refs. [351 EHl EH E] which the reader 
is referred to for more details. The synchronous, time-orthogonal gauge is defined by 
setting g Q0 = —a 2 (r) and = 0, so that the line-element takes the form 

ds 2 = a 2 (r)[-dr 2 + 7ij (x, r)dxW] . (32) 

For our fluid containing irrotational, pressure-less matter plus A, this also implies 
that the fluid four- velocity field is given by m m = (1/a, 0,0,0), so that x represent 
comoving "Lagrangian" coordinates for the fluid element (indeed, the possibility of 
making the synchronous, time-orthogonal gauge choice and comoving gauge choice 
simultaneously is a peculiarity of fluids with vanishing spatial pressure gradients, i.e. 
vanishing acceleration, which holds at any time, i.e. also beyond the linear regime). A 
very efficient way to write down Einstein and continuity equations is to introduce the 
peculiar velocity-gradient tensor [37] 

^ = ^ " -J) = \^ , (33) 

where we have subtracted the isotropic Hubble flow. Here semicolons denote covariant 
differentiation. From the continuity equation T^ v ' . v = 0, we infer the exact solution for 
the density contrast 5 = Sp/p [371 E] 

5(x, r) = (1 + <f (x)) [ 7 (x, r)/ 7o (x)]- 1 / 2 - 1 , (34) 

where 7 = detjij. The subscript "0" denotes the value of quantities evaluated at the 
present time. 
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From Eq. it is evident that in the comoving- synchronous gauge the only 

independent degree of freedom is the spatial metric tensor 7^. The energy constraint 
reads 

tf 2 - Pj&i + me + 11 = m 2 Vt m 5 , (35) 

where 1Z l a is the Ricci tensor associated with the spatial metric 7^ with scalar curvature 
1Z = TZ % i . The momentum constraint reads i?*^ = where bars stand for covariant 
differentiation in the three-space with metric 7^. Finally, one can use the Raychaudhuri 
equation 

& + %■& + & ,^ + ln 2 n rn 5 = , (36) 
j 2 

which is obtained from the energy constraint and the trace of the evolution equation 
+ 2H# i j + &&< j + ~ (tf V k - tf 2 ) S i j + TVj - -KFj = . (37) 

Notice that these equations are exact and describe the fully non-linear evolution of 
cosmological perturbations (up to the time of caustic formation). In order to show 
how the primordial non-Gaussianities appear in in the matter density contrast, we then 
perform a perturbative expansion up to second order in the fluctuations of the metric. 
The spatial metric tensor can be expanded as 

7* = (1 - 2^ (1) - ^ 2) ) * y + xg } + \xf? , (38) 

where xlj an d Xif are traceless tensors and include scalar, vector and tensor 
(gravitational waves) perturbations. As usual we split the density contrast into a linear 
and a second-order part as 5(x, r) = <jW(x, r) + |#( 2 )(x, r). At linear order the growing- 
mode solutions in the comoving-synchronous gauge are given by 

^m (x , T ) = ^ x) + _J_vV(x,x), 

*SW) = Ai)X (1 W), x (1) (x,t) = - — * y(x,T) , (39) 

where </?(x, r) is the growing-mode scalar potential defined in Eq. (II ip . The linear 
density contrast 8^ 1 ' in this gauge is related to ip via the usual Poisson equation, namely 

VV(x,r) = |w 3 (r)fi m (r)(jW(x,r). (40) 

In writing we have eliminated the residual gauge ambiguity of the synchronous 
gauge as in Ref. [6] 0- We have also assumed that linear vector modes are absent, 
since they are not produced in standard mechanisms for the generation of cosmological 
perturbations (as inflation). We have also neglected linear tensor modes, since they play 
a negligible role in LSS formation. 

More in general, at any order n in perturbation theory the scalar potentials ij}^ n ' and x can be 
shifted by arbitrary constant amounts 5i/)q an d &Xo°\ om Y provided Sip^ + (l/6)V 2 5xo^ = 0- 
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By perturbing Eq. 

6 (2) 



r(2) 
°0 



up to second order, we get 
3(^ (2) -4 2) )+^(^- 1 ) [I (V 2 xi 



+ d^x^did^ ~{D + -1 



(i) 



24 1} + 3 V 



2,,M 



v 2 ^ 



(41) 



To compute the metric perturbation ip^ 2 \ we can use the evolution equation, 
Eq. fl37j) . To this aim it proves convenient to write 



(2) 



1 

3 



(2) 



2if>> 



(i) 



v 2 ^ ) v^- 



72_(1) 



1 

12 



v 2 x« 



(42) 



where we have introduced the variable £ which is determined by the equation 



V^ 2 



1 



3e(f2 r 



d^x^dAx'o 



(i) 



(43) 



An approximate solution of this equation can be obtained by making use of the 
usual approximation (see e.g. Ref. (SHUSH]) e(f2 m ) »s 1, which holds true for reasonable 
values of Q m . Under these circumstances we can write 



AD At) 



Dl(r) 
14 



2(V 2 X« 



^x^dAx'o 



(i) 



(44) 



where A is an integration constant which can be determined by the energy constraint 
at the initial time. Notice that it is precisely the sub-leading, Post-Newtonian term 
proportional to the linear growing mode D + which brings all the relevant information 
about primordial and GR-induced non-Gaussianity. It is also important to stress that 
the time dependence of this term, which comes from the homogeneous solution of 
the above equation is exact, while the above approximation only affects the fastest 
growing Newtonian terms, i.e. those proportional to D\. Using this procedure and 
providing the initial data (formally at r — > 0) in terms of the gauge-invariant curvature 
perturbation, which in this gauge reads = — ift^ — (l/6)V 2 xj^ = ~~ 5tp- m /3 and 



.(2) 
'in 



-V> (2) - 
100 



;i/6)v 2 x 



2,,(2) 
in 



+ 



m 2 



50a n i<£> 2 n /9, we finally obtain 
~i -1 



Om) 



-n 



T 



0m y 



0?7I 



3 n 

2 s 'Om 



4 



(2 - a NL )v? in V 2 v2i; 



where we made use of the following property 

3 



UD' 



const. 



(45) 



(46) 



whose validity can be easily proven on the basis of Eq. (|T2|) and of the Friedmann 
equation %' — % 2 + (3/2)'H 2 Q m = 0. One may notice that even the Newtonian part 
(i.e. the one proportional to D 2 + ) of the second-order matter perturbations differs from 
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the Poisson gauge expression. This is a well-known feature that is to be ascribed to 
the different meaning of the mass density when moving from Eulerian to Lagrangian 
coordinates, which only appears at second and higher orders (see, e.g. Ref. [30]). 

As we did in the Poisson-gauge case, we can re-express this result in terms of a 
suitable potential $ in order to introduce a non-linearity parameter /nl with the usual 
meaning. We obtain (see also Ref. [32] ) 



/£ 1 (k ll ka) = | 



/ \ 5 ki 

G NL - 1) - 1 + 



(47) 



2 k 2 

Notice that the non-linearity parameter is different w.r.t. the one in Eq. (l3Tj) . not 
only because they have been computed in two different gauges, but also because in the 
comoving synchronous gauge the Post-Newtonian term giving rise to f^ L can be easily 
written in an exact form where the growing mode D + (t) is factored out, while this is 
not the case for the Poisson gauge. It is also important to stress here that the expression 
for /nl obtained in this gauge is the one to be used to evaluate the effect of NG on the 
Lagrangian bias of dark matter halos, as recently stressed in Ref. [41]. In Ref. [32] it 
was shown that tiny effects in the above equation which come purely from the General 
Relativistic evolution, i.e. the term which survives in the limit cinl = 1, are potentially 
detectable for some planned LSS surveys. 



5. Concluding remarks 

In this paper we have investigated the effect of primordial and GR-induced non- 
Gaussianities on the second-order matter density perturbation in a ACDM cosmology. 
The calculation has been performed in two popular gauges, the Poisson gauge and 
the comoving time-orthogonal one, which are useful for comparison with observations, 
depending on the particular quantity under study. For instance, in evaluating the 
effect of NG on the mass function and Lagrangian bias of dark matter halos, the 
comoving gauge expression is more appropriate, while the Poisson gauge formulae are 
more suitable for gravitational lensing studies. The strongest present limits on /nl come 
from the analysis of the angular bispectrum of WMAP temperature anisotropy data. 
Indeed, Komatsu et al. [12], analyzing the 7-years WMAP data obtain the 95% limits 
-10 < f#f? < 74, and -214 < /^ L uilatcral < 266. The analysis of Planck data both in 
temperature and E-mode polarization is expected to improve the accuracy by almost an 
order of magnitude. A complementary and very powerful information on the amplitude 
and shape of primordial NG will come from the study the galaxy clustering (e.g. 
Ref. [43]) and other LSS datasets, such as weak gravitational lensing (e.g. Refs [441 145] ) 
and redshifted 21cm background anisotropy (e.g. Refs. [461 HZ])- Primordial NG in 
LSS data can be searched for by various techniques: abundance of massive and/or high- 
redshift objects [IE] [49] [50] EE] E2J [531 E3J , abundance of voids [55], higher-order statistics 
such as bispectrum and trispectrum (see, e.g. Ref. [56] ), large-scale clustering of halos, 
thought as rare high peaks of the dark matter distribution [571 E31 EU [58], [59J, [51] , and 
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their cross-correlation with the CMB via the Integrated Sachs- Wolfe effect [50], small- 
scale NG corrections to the matter power-spectrum [591 EE]- For instance, Slosar et 
al. [62], exploiting the scale-dependence of the NG correction to the halo linear bias 
A6 NG oc f^l^k' 2 , obtained the 95% confidence range —29 < f^ 1 < 70, using the 
two-point function of a combination of datasets. The prospects for sensibly narrowing 
these limits with the advent of galaxy surveys sampling regions comparable to the 
Hubble volume are analyzed in Ref. [63]. Indeed, as shown in Ref. [32], there are very 
promising prospects to observe NG signatures down to the limits of the order unity GR 
corrections discussed in this paper. This largely motivates theoretical efforts to obtain 
accurate predictions of these effects. 
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